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<N ; Abstract 

(— I ' In the context of inter-jet energy flow, we present the first quantitative result of the resummation of 

Qh' non-global logarithms at finite Nc- This is achieved by refining Weigert's approach in which the problem 

p ' is reduced to the simulation of associated Langevin dynamics in the space of Wilson lines. We find that, 

O ■ in e'^e" annihilation, the exact result is rather close to the result previously obtained in the large- A'^c 

mean field approximation. However, we observe enormous event-by-event fluctuations in the Langevin 

process which may have significant consequences in hadron collisions. 



> 

o 
en . 

^ . 1 Introduction 

(^ I In certain search channels of the Higgs boson and new particles at the Large Hadron Collider 

m ■ (LHC), it is often desirable to be able to control QCD radiation from tagged jets in order to 

suppress large backgrounds. A prime example is the Higgs boson production in association 

with di-jets. The two competing production mechanisms, the gluon fusion and the vector boson 

/\^ '. fusion processes, have different patterns of soft gluon radiation due to the difference in their 

j^ i color structure. Discriminating these processes quantitatively using some measure of radiation 



is therefore a useful strategy to determine the Higgs couplings II1I2II . 

A related class of observables which are particularly sensitive to soft radiation is the cross 
section with a veto on unwanted jets in the full or partial region of the phase space. This gen- 
erally requires the resummation of logarithms in p^°*°, the threshold transverse momentum of 
vetoed jets. Steady progress in this direction has been made for global observables which in- 
volve all the particles and jets in the final state including those close to the beam axis. The 
state-of-the-art is that one can resum the leading logarithms (LL) (a^ ln^py^*°)", the next- 
to-leading logarithms (NLL) (a^lnp™*")" and even the next-to-next-to-leading logarithms 
(NNLL) (a^inpveto^n ^^ 
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However, in contrast to such progress, there exists a severe limitation in our ability to resum 
non-global logarithms [I5I6I which arise when measurements are restricted to a part of the phase 
space excluding the beam and jet regions. In this case, double-logarithms (a^ln^p™*")" are 
absent due to the lack of the coUinear singularity. The leading contribution is then comprised of 
single-logarithms (as Inp^^")"- which originate from the soft singularity. The problem is that 
these logarithms do not exponentiate, and because of this difficulty their resummation has been 
hitherto done only in the large- A^c limit [I5I6I7I8L In other words, even the leading logarithms 
cannot be fully satisfactorily resummed. This could be a potentially serious drawback in ac- 
tual experiments considering the fact that, strictly speaking, any vetoed cross section at hadron 
colliders is inevitably non-global due to the finite acceptance of detectors. 

In fact, there is a single work by Weigert [9 1 which did discuss the resummation of non-global 
logarithms at finite Nc in a simpler setup of e+e~ annihilation where complicacies from the ini- 
tial state radiation do not arise. His approach is based on an analogy with another, seemingly 
unrelated resummation in QCD, namely, that of the small-x (or 'BFKL') logarithms in Regge 
scattering. In this context, a very similar issue arises as to how one can generalize the equation 
which resums small-x logarithm in the large-iVc limit to one at finite Nc. It turns out that these 
equations bear a striking resemblance to the equation which resums non-global logarithms in 
the large- A^c limit. Since technologies to solve the former problem are well-developed, they 
may be suitably adapted to address the latter problem as well. Somewhat surprisingly, however, 
Weigert's approach has not been pursued for a decade. Part of the reason of this may perhaps be 
that, as we shall point out, there is actually a flaw in his formulation which deters a straightfor- 
ward numerical implementation. In this paper we overcome this difficulty and present the first 
quantitative results of the resummation of non-global logarithms at finite Nc- 

In Section 2, we quickly review the nonlinear evolution equations which resum the small-x 
logarithms in high energy QCD. The subject may seem utterly unfamiliar to the readers whose 
primary interest is jet physics. However, the similarity (or even equivalence) to the resumma- 
tion of non-global logarithms will soon become apparent in Section 3 where we introduce the 
relevant evolution equations. Using this similarity, we discuss how to solve the equation for 
non-global logs at finite N^ in Section 4, and present numerical results in Section 5. Finally, we 
examine the results and conclude in Section 6. 



2 B-JIMWLK equation 

Consider a 'dipole' consisting of a quark and an antiquark located at transverse coordinates 
X and y, respectively. The S'-matrix of the dipole moving in the x^ = 75(2^° — x^) direction 
and scattering off some target in the eikonal approximation is 

{S.y)r = ^JtTiU.Ul))r, (1) 

where Ux is the Wilson line in the fundamental representation 



iU^)ij = Pexp (i r dx-A+{x-,x)t'') , l<i,j <N, 

\ J~oo J ij 



(2) 



and T is the rapidity of the dipole. The expectation value (■ ■ ■ ) is taken in the target wavefunc- 
tion. In deep inelastic scattering (DIS), r is essentially the logarithm of the Bjorken-x variable 

a, , 1 
r = — In - . (3) 

TT X 

The dipole S'-matrix ([T]) (more precisely, 1 — (S)) then measures the total cross section of 
the subprocess YP ~^ QQP — ?• -'f at the corresponding values of x and the photon virtuality 
Qr^l/\x-y\. .,^_^_ 

The B-JIMWLkQ equation 111011111211 resums the small-x logarithms t'' ~ (a^ In l/x)"- 
which arise in the high energy evolution of (S). It reads 

dr{S^y)r = N, j ^M^y{z)[{S^,S,y)^ - {S^y)^^ , (4) 

where the dipole kernel is given by 

M.y{z) = -, ^""'^^^ ., . (5) 

[x — zY[z — yY 

The equation dH) is not closed because the right-hand-side contains the double dipole S'-matrix 
{SS). This non-linearity reflects the gluon saturation effect in the target. The consequence is 
that dH) is actually the first equation of an infinite hierarchy of coupled equations. However, it 
becomes a closed equation — the Batlisky-Kovchegov (BK) equation [I10I13II — if one truncates 
the hierarchy and assumes factorization {SS) — t- {S){S) 

dr{S^y)r = N , j ^M^y{z)[{S^,) r{S.y) r " (S^j;).) • (6) 

While the BK equation ^ can be solved numerically in a straightforward manner, solving 
the B-JIMWLK equation ((H) had been difficult until the ingenious reformulation of the problem 
as random walk [|14,| . To explain this, it is essential to write the equation in the operator form 

C'riSxy/T ^ {-ti Sxy/r ■ {') 

The effective Hamiltonian H takes the form 

72 

H = jd^xdPy^]Cxy{z)Wl (l + iJlUy - m. - UltjyY'wl 

= j d^xd^y'^Kxy{z){l - UiU.r{l - UlUyr^lVl - jd'xalVl , (8) 



^ Acronym for Balitsky, Jalilian-Marian, lancu, McLerran, Weigert, Leonidov and Kovner. 



where f/"'' = (W^'^ is the Wilson line in the adjoint representation and JC is the gluon emission 
kernel 

In the last line of (O, we defined 

< = -^[ ^-^—^tr{T^Um ■ (10) 

J zn [x — zy 

The derivative V acts on the Wilson line in the fundamental representation as 

VlUy = iU^t^5^^\x-y), VlUl = -ieul5^^\x - y) , (11) 

and also on the adjoint Wilson line with {T"-y^ = —ifo.bc 

VlUy = iU^T-5^^^ (x-y), VlUl = -iT^Ul^^^^ {x - y) . (12) 

The latter operation has been used in obtaining the cr-term in ([8]). An important observation 
is that ^ may be viewed as the Fokker-Planck equation treating Wilson lines U as dynami- 
cal variables. As is well-known, there is an associated Langevin equation which describes the 
random walk of these variables. The latter can be simulated numerically on a lattice, and the 
solution of the B-JIMWLK equation has thus been obtained in [15J. We shall later discuss this 
approach in detail. 

Before leaving this section, it is useful to show another form of the Hamiltonian derived in 
1161 



H = -J d'xSy^M.y{z)Vl (l + IJiUy - tjiU. - UlUyf K 

= ^Jd'xd'y^M^yiz) (l + iliUy - UlU, - UlUyY'viV'y . (13) 

This Hamiltonian can be used only when it acts on 'gauge invariant' operators of the form 

tr(f/,f/J), tiiU^UlUMl---), tT{U^Ul)tT{UMl), ■■■ (14) 

and generates the same equations as those obtained from (|8]). Eq. (fT3l) features the dipole kernel 
^ instead of the gluon emission kemel Q. The following relation between the two kemels is 
worth noting 

Ma^yiz) = 2lC^y{z) - /C^3,(Z) - lCyy{z) . (15) 




Fig. 1. Back-to-back jets. 



3 Non-global logs at finite Nc 



We now turn to the resummation of non-global logarithms which is our primary interest. 
Consider, as in the final state of e+e^ annihilation, a pair of jets that is overall color singlet 
and pointing in the direction of the solid angle ria,i3 = {da,i3, 4'a,i3) measured with respect to 
the positive 2;-axis (see Fig. [T]). Let Cm be the region inside a pair of back-to-back cones with 
opening angle 6in which include the jets, and let Cout be its complementary region. We then ask 
what is the probability P{VLa^ Vtp) that the total flow of energy into Cout is less than Eout- In 
the perturbative calculation of P in the regime Q ^ Eout ^ ^qcd where Q is the hard scale, 
logarithms of the form (a^ InQZ-Eout)" appear which have to be resummed. These logarithms 
are non-global because the measurement is done only in Cout- To leading logarithmic accuracy, 
one may identify Eout with the jet veto scale p™*° mentioned in the introduction. 

It has been shown by Banfi, Marchesini and Smye (BMS) that P satisfies the following 
evolution equation \Q\ 



drPa^ = N,J -^Ma^i-f) (Q^n{l)Pa^PyP - Pa^ J , 

where we abbreviated Pa/^ = P{^a,^(s) and the parameter r is defined as 



(16) 



•T Eoul 



llNc-2n 



■ , (fixed coupling) 

^^( in'wAQc^ ) • (ninning coupling) 



(17) 



The equation (fT6l) resums logarithms r" ~ (a^ In EoutY to all orders. The integral kernel 





1 — COS 6al3 


(1- 


- cos6'q,^)(1 — cos6'^,3) 




1 - n« ■ n^ 


(1- 


-n„-n^)(l -n^-np) 




ria -Up 



Map{l) 



is composed of null vectors n'^ = (1, sin 6^ cos 0, sin 9 sin 0, cos 9) = (1, n) proportional to the 
four-vector of hard partons. The 'step function' 

ensures that real gluons are emitted only in C^O A trivial rewriting of (fT&l 

(20) 

where ©out (7) = 1 — Qinil) illuminates the physical meaning of the right-hand-side. The first 
term represents the familiar Sudakov suppression, whereas non-global logarithms are resummed 
by the second term describing the emission of an arbitrary number of gluons into Cj„ which then 
coherently emit the softest gluons into Cout- 

For the purpose of the present work, it is indispensable to note that the equation (fT6l) has 
been derived in the large-A^c limit. Its finite-A^c generalization was discussed by Weigert [|9J 
and the result reads 

dAPa,)r = N^J^M^,{^)(<dUl) [Pa,P,,-^^-'^Pa,'j 

= -2Cf j '^M^p{l)Qout{l){Pap)r 

+^c I ^M^f,i^)Q^nil){Pa^PyP - Pa^) ^ , (21) 

where Cp = -^fw'- However, (|2TI) itself is highly formal in that the meaning of the averaging 
(■ ■ ■ ) can be specified only indirectly, as a proxy of certain complicated functional integrals 
dH. Nevertheless, putting this qualification aside, the striking similarity between Eqs. (fT6] ). (1211 
and Eqs. dU, ^ is unmistakable. In fact, it is possible to establish a rigorous mathematical 
equivalence between the two problems. As shown in [1171181 . a conformal transformation known 
as the stereographic projection 



^ In principle, real gluons can be directly emitted from hard partons into Cout provided their energy is 
less than Eout- However, to leading logarithmic accuracy such contributions may be omitted tli. 



, 1 o, f sin 9 COS (j) sin6'sin0\ 

X = {x\x') = -— ^, J , (22) 

\l + cos 1 + cos J 



exactly maps the respective kernels onto each other 

(fz [x — yY dQ,j 1 — cos 6*0/3 



2tt {x - z^iz - yf An (1 - cose„^)(l - cos^. 



(23) 



Aside from the kinematical constraint factor 0j„, the map (|23] ) dispels any structural difference 
between the two sets of equations |iJ This equivalence probably has a deep geometrical origin 
which goes beyond the perturbative framework. Indeed, such a correspondence persists even in 
the strong coupling limit of A/" = 4 supersymmetric Yang-Mills theory [17J. 

Fortunately, a powerful machinery to solve the JIMWLK problem (H]) is available, and this 
brings hope that the jet problem (|2TI ) can be solved in a similar manner. For this purpose, one 
seeks the operator form of (1211 with a formal identification 

Pa0 ^ j^triU^Ul) . (24) 

Here the Wilson lines Ua,i3 represent eikonalfiet) lines starting from the space-time origin and 
extending to infinity in the direction of fi^ /3[^ Eq. (|2TI) can then be written as a Fokker-Planck 
equation in group space 

dr{Po.p)r = -{HP^p), (25) 

where BH 

= ]^ j d^^dnp'^Map^i) (l + UiUp - Qin{l){uiU, + U\Up)Y\l% . (26) 
In ([261), the derivative V is defined by {diVt - Vt') = 6{cos 6 - cos 6'')5(0 - 0')) 

VlUp = iU^e5{n^ - Up) , VlUl = -it'^Uidin^ - np) , (27) 

and similarly to (fT2]) in the adjoint case (f^ — )■ T°-). 



^ Note that if one sets 0j„ — )• 1, (|2TI ) becomes identical to ^ under the map (|23] ). 
^ Given the probabilistic nature of P, one should more properly consider U as the product of a Wilson 
line in the amplitude and that in the complex-conjugate amplitude [9]. In practice, however, this does 
not matter in the equivalent Langevin approach. 



4 Equivalent Langevin dynamics 

The Fokker-Planck equations (|7]) and (|25] ) can be solved by making use of an equivalent 
Langevin formulation. To illustrate the idea, consider the following Fokker-Planck equation for 
some probability distribution Pt{x) of dynamical variables {x""} 



dPJx) 1 



dr 



-da (r\x)dbPr{x)) (28) 

\x''\x)dadb + (y''d^P,{x) (29) 



■da 



h,(r\x)Pr{x))-a^Pr^ 



X] 



(30) 



where x"^ = x^° and a"- = ^d^x''^- We assume that x is factorized in the form x"'' = S°-'^£'^^. 
The equivalent Langevin equation is then 

^x'^{r) = a''{x)+8''%x)Ur). (31) 

dr 

where ^ is the Gaussian white noise characterized by the correlator 

{Ur)Ur')) = 5abS{r - r') . (32) 

The distribution can then be obtained by averaging over an ensemble {a;"(r)} of random walk 
trajectories 

Pr{x) = {5{x - x{t))) . (33) 

To be more precise, the equation (ISTT ) makes sense only in a r-discretized form. There is 
a well-known ambiguity in how we discretize the equation, the so-called Ito-Stratonovich 
dilemma. The appropriate choice corresponding to (|30l l is the Ito scheme 

x-{t + e) = x\t) + ea\x) + ^£''\x{T))Ur) , (34) 

where e is the time step and the argument of 8 is evaluated at the previous time r, respecting 
causality lfT9l . In (|34|) . we have rescaled the noise as y/e^ — > ^ in order to make explicit the 
fact that the typical variation Ax is 0{-/e) in a random walk. With this normalization, the noise 
correlator in discrete time reads 

{Ur)Ur')) = SabSrr' • (35) 

The operator form of the B-JIMWLK equation ^ has precisely the structure (|30] ) with a 
factorized kernel ([8]). It can thus be described by the Langevin dynamics (|34l) . with the SU(3) 
matrices Ux,y playing the role of {x'^} [fT4l . One can simulate the random walk in group space 



on a transverse lattice, and this is how the solution to B-JIMWLK equation has been originally 
obtained lfT5ll. 

In [P], Weigert suggested to follow the same strategy in solving (|25] ). A comparison of (|26l ) 
and (|30)) implies that a = 0. The kemel of (|26]) 



= Qout{l) + QoutilWiilp + Q^n{l){l ' UiU^){l - UlU p) , (36) 

can be written as a sum of three factorized terms. [Note that (0m)^ = ©in-] Exploiting the 
factorized form of M. as seen in the last line of (fT8l) . Weigert deduced an analog of the stochastic 
term 8"'''^c in (l3T| ) by introducing three independent noises ^^^\ (/ = 1, 2, 3) 



£''%c - 1 dn,^^^ {eui) {s^^l + (uir^l) + 0.n(7)(i - uiu.r^^l} , 

(37) 



Ct ' ' ^'y 



which generates a random walk U — )► Ue^^"^""^" in group space. However, the fact that A^ is 
factorized in four-vector space means that the noises must have the correlator 

(^(/)M^(^)-)^5^^5/J^M^^ (38) 

which is negative for the spatial components fx,^ = 1,2,3, hence they cannot be simulated in 
practice. 

The resolution of this problem again comes from the correspondence with the B-JIMWLK 
evolution. Firstly, the identity (|23T ) clearly shows the two-dimensional nature of the dipole ker- 
nel, so introducing four-vectors is an excess. We then notice the close similarity between (|26l) 
and (fT3l) . the latter being equivalent to dS]). This suggests that we can rewrite the effective Hamil- 
tonian in a form analogous to ([8]) also for the jet problem. We thus look for a kemel /C satisfying 
(cf. ([13) 

Ma^ij) = 2/C„;3(7) - /C„„(7) - /C;3/3(7) • (39) 

The solution is found to be (cf. dH)) 

cos 9ay + COS 9^^ - COS 9ai3 - 1 _ (n^ - n^) ■ (n-y - n^) 
2(1 — cos6'a7)(l — cos6'^,3) 2(1 — n„ • n.^) (1 — u-y • U/j) ' 

which has a factorized structure on the unit sphere (|n| = 1) embedded in three spatial dimen- 
sions. Reversing the argument in [16] which led ^ to ([TJt . we arrive at an effective Hamiltonian 
equivalent to (|26l) 



= I dn^dn^^jc^f^i^) (i + uiUf, - eui){ulu-y + u\Up)Y\lv''p 

-jdn^alVl, (41) 

where 

< = -z / f^ ,, ®^"^^^ M T'^iJiU,) . (42) 

J 47r (1 — riQ ■ n^) 

Moreover, differently from (|36^ . we write the expression in the brackets as a sum of two factor- 
ized terms 



= (1 - QUl)UiU,){l - Q^nU\Up) + QoutilWitjp . (43) 

This reduces the number of independent noises from three to two ^^^^ (/ = 1, 2). They are 
characterized by the following correlator (in discretized 'time' r, cf. (|35] )) 

(ei?(r)4f (r )) = 5..'5(cos^, - cos^;3)5(0„ - (^p)5"5ab5^' , (44) 

where k^l = 1, 2, 3 are the spatial indices. The problem of negative metric has been circum- 
vented. We can now write down the associated Langevin evolution for Ua with Via G Cj„ (cf. 

U^ {r + e) = f/. (r) exp [i^ [v~s J dQ, (£«- • |« + ^i^^^ ■ ^g)) + e<) } , (45) 



where 



i^T)' = ^r^V^fi - ^-(^)uiu,r , (46) 

V47r i — no, ■ n^ 



(4t^)' = -^ i''" ""f e..,(7)(?7t)- . (47) 

' v47r 1 - tIq, ■ n^ 



All the matrices on the right-hand-side of (1481 ) are evaluated at r according to the ltd scheme. 
Expanding the exponential up to 0(e), we get 



u^T + e)^ u^{t) + tViu^{T)e J dn, (£:«^^ ■ |« + 4^)'^^ ■ |( 

+^f/a(r) |^r< - i (f^Jdn, (£«- ■ |W + £:g- ■ 42,)))'! . (48) 



10 



The second line of (|48l) can be simplified as follows. In the cr-term, we use the identities U"-''t^ 
WrU and rf/r = |(trf/) -^Uto obtain 



[/„t"tr(T"t/t?7^) = -f/„[t^t^](f/t)-'^(t/^; 



db 



= itr(f/,f/t) U^ - itr(f/^f/t) f/,f/tf/„ . (49) 

To the accuracy of 0{e), the terms quadratic in noise ^^ may be replaced by their expectation 
values using (|44)) . After these manipulations, (1481 ) takes the form 

\k 



+^ / T^l ^ f-2C^t/a + e,„(7) ftr(f/,t/t) f/^ - -^f/.)) . (50) 

J 47r 1 — no, • n^ V V ' ISc / / 

Note that there is no singularity at fi^ = Qq, G Cin- By computing the difference 

i-tr(t/,(r + e)Ul{T + 5)) - i-tr([/Jr)f/t(r)) , (51) 

to 0{e) and using (|44|) . (|39| ) and the relation /CaQ,(7) = —1/(1 — n^ ■ n^), one can recover (1211 ) 
after the identification ([24)) . 

For the sake of numerical simulations, it is more economical to express the evolution (|50] ) in 
a left-right symmetric fomo 



t/„(r + e) = e*^"'f/,(r)e^^", (52) 



where 



^^ = \l]^ J d^^ ^C- n'^.n {-^^MU.t^Ui^!,'^' + QUl) n^^) , (53) 



"a '"7 



s f r^ (nry — n^ ''^ 



^"-^Tj"''^T^7d^''^^^'' <=^' 






In this representation only the terms proportional to noise are kept in the exponential. It is easy 
to check that (|52l) and (l50l) are equivalent to 0{e) under the identification ^^ ^ (^^). Eqs. (l44l) 
and (|52l) will serve as the starting point of our numerical simulation. 



^ See Ref. |[20l for a similar rewriting of the JIMWLK evolution. As noted in this paper, the right- 
multiplication rule such as (1451 ) is related to the choice of the Hamiltonian ^ or (l26l ) being expressed 
in terms of 'right-derivatives' V^[/ ~ Ut"-. The cr-term can be eliminated in the process of converting 
one of the right-derivatives into a left-derivative VJC/ ~ t"-U. 
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5 Numerical simulation 

We simulate the random walk (|52] ) of Wilson lines living on the unit sphere by discretizing 
the coordinates 1 > cos 9 > —1 and 2n > (f) > with lattice spacings Oc and a^, respectively. 
The SU(3) matrices Ua are defined only at the grid points belonging to Ci„, whereas the noises 
C,i^^ are defined at all grid points [^ The initial condition at r = is simply given by 

Ua = I , (unit matrix) (55) 

for all Via E Cin, or equivalently, P^^ = 1 for all pairs {Via, ^/s) corresponding to no radiation 
before evolution]^ We then update {Ua} after each time step e according to the formula (|52] ) 
with noises ^^^^ (/ = 1,2) randomly generated from the Gaussian distribution 

n yW^^P (-"-^^'.'^%T) , (ei?4f ) = -^S^'6.,6-6a, . (56) 

In order to ensure that f/'s remain unitary during the evolution, we need to evaluate the expo- 
nential of matrices e*"^^/^ accurately (although the equation (|52] ) makes sense only to 0{e)). In 

l+2^ + """ + ;^(2^) ) with m,r2 
large enough. On top of this, we perform the 'reunitarization' of t/'s using polar decomposition 
method after every 100 steps of evolution. The effect of this latter operation is actually very 
small due to our accurate evaluation of e*"^^/^. 

The above procedure is repeated a desired number of times A^ = r/e, and at the end of this 
random walk trajectory we compute the trace 



tr(UJT)U'' 



tTiUaiT)U}ir)) . (57) 



c 



We then average it over many such trajectories and identify the result with {Pai3)T- In practice, 
we choose e = 10^^ and average over 500 independent trajectories. 

Fig. |2] shows the evolution of P^ (^ajj^/?) for cos 9a = ^ and cos 9j3 = —1, corresponding to 
back-to-back jets in the beam direction Lj The opening angle of the cones (see Fig. [I]) is fixed to 
cos 9in = ■^. The simulation was done on a lattice with 80 grid points in the cos 9 direction, and 
40 grid points in the direction. The exact Nc = 3 solution to (|2T1) (solid red line) is compared 



^ At cos 6a = ±1, Ua and ^a are independent of (j), as they should. In the case of Ua, this is guaranteed 

by our initial condition ( [55] ) and the structure of the evolution dSOl l which preserves this property. 

^ In the IIMWLK case, the initial condition is the value of the 5-matrix ([T) at small, but not too small 

value of X. This is model dependent and its initial sampling is non-trivial ifTSll . 

® We actually plot the real part of the average {tiUaUl). tiUaUl is in general complex-valued for each 

trajectory, but we have checked that the imaginary part of the average is consistent with zero within 

errors. 
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Fig. 2. Solid line (red): exact Nc = 3 solution to (|2TI ). The band indicates the standard error. Dashed line 
(blue): Nc = 3, mean-field solution to (1581 ). Dotted line (green): solution to the BMS equation (fT6l ) from 
Elil . Dash-dotted line (yellow): result with only the Sudakov teiTn. 



with the solution of the large- A^c BMS equation (1161) (dotted green line) previously obtained in 
[|21llJ^| and also with the solution of the 'mean field approximation' to (|2TI) (dashed blue line) 



dr{Pap)r = -2C p J -^Mo.p{j)eoutil) (Pa^) r 



+iVc 



An 



■Map{-r)Qin{l){{Pa^)r{P^p)r " (^a/3)r) , 



(58) 



which differs from the BMS equation only by the coefficient of the Sudakov term Nc = 3 -^ 
2Cf = 8/3. The latter serves as an indicator of the quality of the mean field approximation 
{PP) — )■ (-P) (P) . For the sake of reference, we also plot the solution obtained by keeping only 
the Sudakov term (first term on the right-hand-side) in (|58l l (dash-dotted yellow line). 



^ Note that the definition of r in |[2ll differs from dTT]) by a factor of Nc 
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6 Discussion 

A comparison of the solid (red) and dashed (blue) lines in Fig.[2]shows that the exact solution 
is rather close to the mean field solution with the finite- A^c corrected Sudakov term for all 
values of r explored in this work. This may come as no surprise to those who are acquainted 
with the solution of the JIMWLK equation which agrees well with the BK solution II15I22II . 
Nevertheless, we find the present result quite intriguing because we actually observed enormous 
trajectory-by-trajectory fluctuations. In Fig. |3]we show 50 (out of 500) individual random walk 
trajectories used in the computation of the average. It turns out that the fluctuations are so large 



that the standard deviation 5P = J (PP) — {P){P) is of the order of (P) itself for not-so- 
small values of r. Actually, this is the reason why we needed to run (9(100) trajectories to 
obtain a reasonably stable result. Such large fluctuations have not been seen in the previous 
simulation of the JIMWLK equation where 'already one trajectory gives a good estimate of the 
final result' [15J. 

In our opinion, the crucial difference between the two problems which has resulted in such 
different behaviors of fluctuations is the initial condition. In the jet problem, the initial condition 
{P)r=o = 1 means that there are no partons besides the qq pair. In the parlance of saturation 
physics, the system is initially very 'dilute'. The evolution then produces soft gluons which 
multiply exponentially in r according to the BFKL formula [I8I17II . It has been demonstrated 
that these gluons have very strong number fluctuations [|23l and spatial correlations [24 25 18 1. 
We thus find it natural to attribute the observed large values of 5P to such fluctuations and 
correlations. On the other hand, when solving the JIMWLK equation, one often uses 'dense' or 
'classically saturated' initial conditions; {Sxy)T=o ^ 1 for las — ■y| larger than some value. Since 
there are many uncorrected gluons in the system from the beginning, there is little room left for 
pure-BFKL evolution, i.e., it is suppressed by the saturation effect. Accordingly, fluctuations 
and correlations can develop only weakly, and this is consistent with what has been found in the 
previous simulations. 

In conclusion, we have demonstrated for the first time the resummation of non-global log- 
arithms at finite N^ to leading logarithmic accuracy. Our study shows that, at least in e+e^ 
annihilation and for phenomenologically interesting values ofr < 0.2 ~ 0.3, event-averaged 
non-global observables may be reliably computed in the mean field approximation by solving 
the (modified) BMS equation (f58l) . or equivalently, by Monte Carlo simulations [5^. However, 
the observed large fluctuations imply that the situation may be drastically different for hadron- 
hadron collisions (cf. [26J). Since four partons are involved in hard scattering, one has to deal 
with multiple products of Wilson lines such as tr(f/f/t)tr(f/f/t) and tr (UWUW) (cf. [27 1). 
Moreover, if there are gluons in the initial and final states, each of them picks up an adjoint Wil- 
son line U which further increases the number of (fundamental) Wilson lines. [Roughly, U acts 
like the square of U.] The proper treatment of fluctuations laid out in this paper is potentially 
very important for such observables. We leave this problem for future work. 
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CD 

DC 



large-N 

finite-Ng (mean field approx.) 

finite-Ng (exact) 




0.3 0.4 0.5 



Fig. 3. Thin lines (pink): 50 random walk trajectories used in the computation of the average (thick red 
line). 
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